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ABSTRACT 

We propose that a significant fraction of the ultracompact H n regions found in massive 
star-forming clouds are the result of the interaction of the wind and ionizing radiation 
from a young massive star with the clumpy molecular cloud gas in its neighbourhood. 
Distributed mass loading in the flow allows the compact nebulae to be long-lived. In 
this paper, we discuss a particularly simple case, in which the flow in the Hii region 
is everywhere supersonic. The line profiles predicted for this model are highly charac¬ 
teristic, for the case of uniform mass loading. We discuss briefly other observational 
diagnostics of these models. 
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1 INTRODUCTION 

The formation of massive stars is an important, but vexed, 
question, largely because such stars disrupt their natal en¬ 
vironment by their powerful winds and ultraviolet radiation 
fields. The ultracompact Hu regions (UCHuR), which are 
found deep within molecular clouds, may provide significant 
information on the properties and environment of the very 
youngest massive stars. UCHuR. are highly obscured at opti¬ 
cal wavelengths and most information on them derives from 
infrared and high frequency radio data. A comprehensive 
overview of their properties is given by Churchwell (1990). 
They have high emission measures (( n^)L ^ 10 7 cm” 6 pc) 
and small scale sizes (L ^ 10 17 cm) and thus high r.m.s. 
electron densities ((rig) 1 / 2 <; 10 4 cm” 3 ). The ionizing photon 
supply rates necessary to maintain the observed ionization 
are in the approximate range 10 44 10 49 s _1 , corresponding 
roughly to ZAMS spectral types B2-05 for single star ion¬ 
ization. There is, however, considerable evidence from com¬ 
parison of IR and radio data, that groups of stars may be 
embedded within a given UCHuR and also that the dust 
content of the region may absorb nearly all the UV photons 
generated by the stars (Kurtz, Churchwell & Wood 1994). 
Consequently, the assignment of a unique spectral type to 
the exciting source is an extremely uncertain procedure. 

UCHuR have a wide range of observed morpholo¬ 
gies. Churchwell (1990) defines 5 main morphological types: 
cometary (~ 20 per cent), core-halo (~ 16 per cent), shell 
(~ 4 per cent), irregular or multiply peaked (~ 17 per cent) 
and spherical and unresolved (~ 43 per cent). Most theo¬ 
retical attention has been given to cometary regions (e.g. 
Van Buren & Mac Low 1992). They can be modelled fairly 


convincingly as the steady-state partially ionized structures 
behind bow shocks driven into molecular cloud material by 
the winds of moving stars. 

Although Mac Low et al. (1991) suggest that morpholo¬ 
gies other than cometary (in particular core-halo) can be 
explained as cometary structures viewed along the axis of 
symmetry, it seems likely that other models need consider¬ 
ation. For example, Kurtz, Churchwell & Wood (1994) give 
evidence for a size-density relationship for spherical or unre¬ 
solved regions. This is neither present nor would be expected 
for cometary UCHuR. Dyson (1994) suggested that a nat¬ 
ural explanation for at least some UCHuR arises from the 
interaction between an early type star and a very clumpy, as 
opposed to relatively homogeneous, molecular cloud. Molec¬ 
ular clouds are well established to be clumpy on length scales 
down to the limits of observational resolution. Clumps can 
act as localised reservoirs of gas which can be injected into 
the surroundings by photoionization and/or hydrodynamic 
ablation. The continuous injection of material into a Hii re¬ 
gion can lead to a quasi steady-state configuration which is 
bounded by a recombination front (RF) as opposed to the 
ionization front (IF) which bounds standard Hii regions. 
This particular configuration avoids the expansion problems 
associated with a ‘classical’ Hii region model for UCHuR, 
which would expand on a timescale far shorter than the 
characteristic 10 5 ” 6 yr lifetime of UCHuR estimated from 
the statistics of UCHuR association with early type stars 
(Churchwell 1990). Hn regions bounded by recombination 
fronts can be held in equilibrium by gas pressure or recoil 
pressure depending on whether the flow exits sub- or super¬ 
sonically. 

We discuss here a particularly simple model where a 
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star with both a UV radiation field and a powerful hyper¬ 
sonic wind interacts with clumpy material and show that 
it may explain the structure of shell-like UCHuR. Further 
papers in this series will investigate in detail the full range 
of models outlined by Dyson (1994). A future key test of all 
these models will be the observational determination of the 
kinematics of ionized gas, molecular gas and neutral hydro¬ 
gen associated with such UCHuR. 

Non-steady models with photoionization-induced mass 
injection have also been suggested by Lizano et al. (1995). 
One major difference between their models and those given 
here and by Dyson (1994) is that Lizano et al. assume a very 
specific model for the mass injection process which leads to 
a close coupling between the spatial and temporal behaviour 
of the mass injection and the radiation field. We treat mass 
injection as a free parameter and assume here that it is con¬ 
stant in space and time. Dyson (1994) has discussed some 
of the factors involved in the mass injection process and em¬ 
phasised that both hydrodynamics and photoinjection may 
play roles. It is clear that both approaches need investiga¬ 
tion. 


2 UCHuR AS WIND DRIVEN FLOWS IN 
CLUMPY MOLECULAR CLOUDS 


In this initial paper, we assume that the hypersonic wind 
from a single star blows into a clumpy molecular cloud and 
that mass injection into the flow occurs at a constant volume 
rate q near to the star. Mass injected into a fast wind will 
slow the wind down. The frictional energy dissipated in this 
process may be radiated away by, for example, enhanced ra¬ 
diative losses in boundary layers (Hartquist & Dyson 1993). 
We assume that this occurs here and that the flow is pho- 
toionized by the stellar radiation field and remains isother¬ 
mal at the usual photoionization equilibrium temperature 
T ~ 10 4 K. A more detailed discussion of the isothermality 
of mass loaded winds is given elsewhere (Williams, Hartquist 
& Dyson 1995). 

We assume that mass loading occurs only in the ion¬ 
ized region and that the flow always remains supersonic. We 
neglect the possibility of global shocks in the flow produced 
by mass injection. This will be a valid assumption provided 
that the Mach number in the ionized region is not predicted 
to fall below about 2 (Williams et al. 1995). The recombi¬ 
nation front which then occurs has supersonic inflow and 
outflow (cf. Newman & Axford 1968). A discussion of the 
consequences of global shocks and subsequent transonic flow 
for RF structures is given elsewhere (Williams & Dyson, in 
preparation). 

For steady supersonic flow with uniform mass loading, 
the mass and momentum conservation equations are respec¬ 
tively 


pu 





d 


4nr 2 pu 2 = M*u* = /i*, 


( 2 ) 


where u and p are respectively the flow velocity and density, 
r is the radial coordinate and M* and v* the stellar wind 
mass loss rate and velocity. We have assumed that the mass 


flux in the flow is dominated by gas ablated from the clumps. 
Then from equations ([!]) and (^), 
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Clearly, radiation such as bremsstrahlung and recombina¬ 
tion lines, which have emissivities proportional to p , will 
be dominated by the outer regions of the flow. 

The RF occurs at a radius, Rr, where the photon out¬ 
put rate balances the recombination rate. We assume an 
effective stellar output of S* photons s ~ 4 in the Lyman con¬ 
tinuum (which may be only a small fraction of that actu¬ 
ally generated by the star because of absorption by dust). 
We assume (and justify later) that if mass injection is a re¬ 
sult of photoionization, it needs only a negligibly small frac¬ 
tion of the available UV. If the observed nebula is ionization 
bounded, then it must be powered by a Lyman continuum 
flux of 


rR R 
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4ivr' nfadr 1 , 
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where n (= p/fh) follows from equation (|s|), m is the mean 
mass per nucleon and @2 (= 2 xlCU 13 cm 3 s _1 ) is the hydro¬ 
gen recombination coefficient. Thus from equation (|jj) 
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We treat q as a free parameter since it is presumably deter¬ 
mined by factors such as the density of mass loading centres 
and the mode of ablation (Section 1). 

Equation (Q) implies that we can neglect any apprecia¬ 
ble contributions to the radio emission from the mass loading 
centres themselves (e.g. from dense ionized clump surfaces), 
and so the measured r.m.s. electron density is that of the 
flow and is given by 
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The Mach number just before the RF is 



3/i* 

4-KqR^c\ 


( 8 ) 


Inserting characteristic values S 48 = S*/10 48 s -1 ; ( 12 s = 
/i*/10 28 gems -2 ; R 17 = f?.R/T0 17 cm; Ci = lOkms -1 and 
fh = 2xl0 -24 g, equations (j|7j) and (||) become 

(n 2 ) 1/2 = 3.5xl0 4 S^ 2 R 1 ' 7 3 / 2 cm -3 (9) 

Mr = 0.78S 4 - 1/4 R- 1/4 A28 2 - (10) 


The mass injection rate required is 

q = 2.4xlO“ 30 /i 2 8Rr7 3 -Mr 1 gcm ~ 3 s -1 . (11) 


The Mach number at the recombination front from 
equation ( 10 ) is smaller than required, for the values of S 48 
and R 17 derived from observation, unless the momentum in¬ 
put (for the observed S 48 ) is considerably larger than that 
found for field OB stars. Dust significantly reduces the frac¬ 
tion of S* which ionizes the nebula (and thus increases the 
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effective however, in the sample of Kurtz et al., the 

fraction of the stellar luminosity absorbed by dust is less for 
the earlier type stars which would have stronger winds. The 
presence of multiple stars within a single UCHllR may also 
have an effect. 

It is certain that the flow will have significant density 
inhomogeneities. These can decrease the radius of the re¬ 
combination front for given mass fluxes and ionizing lumi¬ 
nosity. We assume that the inhomogeneities have a large 
density contrast, and that the dense gas dominates both the 
emission and the mass flux compared to the more tenuous 
surrounding gas. In this case, the Mach number of the re¬ 
combination front [equation (10)] scales as Mr oc 
where / is the volume filling factor of the dense regions; the 
mean density of the flow at any radius is / 1//2 times smaller 
than that derived (for observed S'* and Rb) by assuming a 
smooth flow. High resolution images of UCHllR do indeed 
show a strongly clumped structure (e.g. Kurtz et al. 1994). 
In the present paper we assume that the dense gas has a large 
(cs 1) covering factor as seen from the central star, and that 
the dense gas is well coupled to the global flow. Exploration 
of these assumptions demands a treatment of intermediate 
scale structure in the flow (Dyson, Hartquist & Williams, 
in preparation). We have also assumed here that the filling 
factor is independent of radius. 

To see whether the required mass injection is plau¬ 
sible and also compatible with the assumptions in the 
model we treat the very special case where the mass in¬ 
jection comes from the surfaces of pressure confined self- 
gravitating isothermal clumps (Dyson 1968; Kahn 1969). 
Kahn (1969) showed that the mass loss rate from such a 
photoionized globule had the remarkable property of being 
dependent only on the sound speed c n in the neutral gas 
of the globule. We adopt his mass loss rate from a globule 
m g ~ 10 2o c|gs _1 , where the sound speed in the cloud is 
c n = 0.3c3kms _1 to allow for its being cold and molecu¬ 
lar. The number of clumps N c required within the UCHllR 
is thus N c ~ {4n/3)Rlq/m s ~ lSOR 1 ^ 4 f 1 / 4 S ^ 4 p 1 ^ 2 % 4 
which appears plausible. The fraction, <j>, of photons used up 
in the photoinjection process itself is <j> = (47rI?R(j/3mS*) ~ 
6x 10 “ 3 R^ 4 / 1/4 S 4 8 3/4 A 2 g 2 , which is consistent with the as¬ 
sumptions made. 

Mass loss driven by ionization from gravitationally con¬ 
fined clumps is just one of a wide range of possibilities. Hy- 
drodynamical ablation may play an important role in re¬ 
leasing gas from the clumps, for instance, as may magnetic 
fields in maintaining inhomogeneities within a global flow. 
There must be several tens of mass loading clumps within 
the UCHllR, if our continuum approximation is to be reason¬ 
ably valid. There is strong observational evidence for such 
clumping in the emitting gas, ten or more individual emis¬ 
sion peaks being visible in the better-resolved images in the 
survey by Kurtz, Churchwell & Wood (1994). 

This characteristic number of clumps implies a space 
density of ~ 10 6 pc~ 3 . The clumps may be related to the 
partially ionized globules (PIGS) and protoplanetary disks 
around low mass stars (proplyds, O’Dell & Wen 1994) ob¬ 
served distributed in Orion nebula. The peak number den¬ 
sity of stars in the low mass stellar cluster around the 
Trapezium OB stars is at least 5xl0 4 pc -3 (McCaugh- 
rean & Stauffer 1995). While the observed velocity struc¬ 
ture of these sources (Massey & Meaburn 1995) may well 



Figure 1 . Emission measure as a function of z (fractional offset 
from the centre of the UCHllR) for Mr = 2.2, normalized to unit 
total emission. 

be similar to that of the mass-loading clumps in UCHllR, 
the Orion nebula is a far older Hu region than those dis¬ 
cussed here. It is likely that in the near vicinity of a very 
young massive star, there will also be many, more transient, 
local density enhancements which have not quite passed 
the threshold for star formation. In 10 5 yr, between 1 and 
10 Mq of gas will pass through a near-sonic recombination 
front, suggesting a mean density of clump gas in the region 
(3xl0 5 — 3xl0 6 cm -3 ), comparable to the mean mass den¬ 
sity of the proplyd stars in Orion. 

It is also interesting to compare the Keplerian veloc¬ 
ity at the edge of an UCHllR (around lkms -1 ) and the 
relative velocity necessary to transit the UCHllR in its life¬ 
time (around 0.3kms _1 ) to the mean velocity dispersion of 
the Trapezium stars (around 3kms _1 ). A significant flux of 
new clumps will enter the UCHllR during its lifetime, either 
through random motions or systematic infall; the density of 
clumps within the region may also be significantly enhanced 
by gravitational focussing effects. 


3 A SIMPLE MODEL UCHuR 

In order to satisfy the requirement that no global shocks 
occur in the UCHuR, we require Mr ^ 2 (Section 2), 
i.e. {fS asRit/U% s) ^ 0.03. We therefore choose as an il¬ 
lustrative set of parameters / = 0.02, S 48 = 1, R 17 = 1, 
p 28 = 1. We then find the observed r.m.s. density, (n 2 ) 1 / 2 — 
3.5xl0 4 cm -3 (equivalent to an r.m.s. density in the 
smoothed flow of 5xl0 3 cm -3 ), q ~ 1.2x 10~ 30 gcm -3 s _1 , 
<j> = 2.3x 10~ 3 and Mr = 2.2. 

In Figure |l], we show the emission measure, EM = 
f n 2 dl, as a function of offset z from the star (where z is 
given in units of Rr). The total emission is normalized so 
that 

f EM(z)27rz dz = 1. (12) 

Jo 

The rapid rise of density close to the recombination front, 
cf. equation ([|), leads to an edge brightened structure in 
both line flux and in the far stronger free-free continuum. 
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(a) (b) 




Figure 2. Line profiles at (a) 2 = 0.1 and (b) 2 = 0.7, normalized to unit total emission across the nebula (when integrated over the 
area of the nebula in units of Hr and projected velocity in units of c;), for Mr = 2.2. 



Figure 3. The structure of an isothermal model UCH 11 R, for 
parameters given in the text. The solid curve shows the Mach 
number of the flow, the dashed curve the flow density (in units of 
100cm -3 ), and the dotted curve the ionized fraction, x. 


The peak intensity is more than twice that at the centre of 
the of the region, for Mr = 2 .2. 

In Figure |s], we show line profiles for an optically thin 
recombination line (/ oc n 2 ) at 2 = 0.1 and 0.7. For flows 
with Mach numbers m <; 2 at the recombination front, the 
predicted line profiles are symmetric but noticeably double- 
peaked, particularly in the weaker emission lines close to the 
centre of the nebula. For smaller Mach numbers at the RF, 
the line profiles will be singly peaked, although noticeably 
broadened and variable across an individual UCH 11 R. Con¬ 
tinuum optical depths may lead to a systematic blueshift of 
the emission in lower frequency recombination lines. Clearly, 
observations of radio recombination line profiles is a key test 
of the model presented here. 

In Figure we show the structure of the UCH 11 R calcu¬ 
lated for the simple model parameters given above, assum¬ 
ing that the gas is everywhere isothermal, at T = 8000 K (cf. 
Arthur, Dyson & Hartquist 1994). This is a fairly reasonable 


assumption within the ionized nebula. Outside the recom¬ 
bination front, the gas would be expected to cool adiabati- 
cally; however, since the gas is supersonic here, the cooling 
will have little dynamical effect. Beyond the ionized region, 
the Mach number plotted should be interpreted as the flow 
velocity scaled to the isothermal sound speed. Note that 
the neutral hydrogen emerging from the RF has a veloc¬ 
ity of 45 kms . Clearly, the observation of neutral material 
with high radial velocities close to the projected centre of 
an UCH 11 R would again constitute a crucial test of such 
models. 


4 DISCUSSION 

In this paper, we have discussed one of the mechanisms 
for the formation of ultracompact Hii regions introduced 
by Dyson (1994). Winds around young massive stars mass- 
loaded by the remaining shreds of the molecular cloud from 
which the star initially formed can naturally explain many 
of the properties of spherical UCH 11 R. 

For parameters characteristic of spherical UCH 11 R, the 
mass-loaded stellar wind will remain supersonic at all radii, 
so long as the radiating gas has a filling factor / ^ 0.02. 
This filling factor is close to the preferred range 0.03 — 0.10 
quoted by Afllerbach et al. (1994) from a non-LTE analy¬ 
sis of recombination line ratios for the (cometary) UCH 11 R. 
G29.96—0.02, although limb brightening on a finer scale 
than the observational resolution would bias their results 
towards small values. 

We have shown predicted line profiles and distribution 
of emission measure across and UCH 11 R. The double-peaks 
of line profiles, particularly close to the centre of the neb¬ 
ula, are a strong observational diagnostic, so long as suffi¬ 
cient spatial and spectral resolution can be obtained. Ob¬ 
servations of sufficient quality are now beginning to become 
available for analysis. 

In future papers, we will extend the work presented here 
to treat centre-bright nebulae (most naturally understood 
as flows which are subsonic for a substantial range of radii). 
We will also include more details of the radiation transfer 
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process (such as continuum optical depths and stimulated 
emission), of the wind-clump interfaces and of the distribu¬ 
tion and dynamical effects of dust. Preliminary models in¬ 
corporating more realistic heating and cooling rates predict 
temperatures appreciably greater than 10 4 K in the inner 
regions. This will result in edge brightening even more pro¬ 
nounced than that shown above (Williams, Dyson & Red¬ 
man, in preparation). 

Radio recombination lines from helium and carbon, 
while fainter than those of hydrogen will give better veloc¬ 
ity resolution; the differences in ionization potential between 
these atoms means that comparison of their line profiles will 
yield important additional information on the radial struc¬ 
ture of the nebula. 
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